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Abstract. The scheme developed by Hartle for describing slowly rotating bodies 
in 1967 was applied to the simple model of constant density by Chandrasekhar and 
Miller in 1974. The pivotal equation one has to solve turns out to be one of Heun's 
equations. After a brief discussion of this equation and the chances of finding a 
closed form solution, a quickly converging series solution of it is presented. A 
comparison with numerical solutions of the full Einstein equations allows one to 
truncate the series at an order appropriate to the slow rotation approximation. 
The truncated solution is then used to provide explicit expressions for the metric. 
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1. Introduction 

Until now, no one has succeeded in finding an analytic solution to Einstein's equations 
that describes an isolated, rotating, three dimensional perfect fluid. If such a solution 
is to be found, then most likely for a homogeneous body in equilibrium. After all, 
three limiting cases are known analytically for this equation of state: the non-rotating 
limit (inner and outer Schwarzschild solution), the portion of the Newtonian limit 
made up of the Maclaurin spheroids and the disc limit (the relativistic disc of dust, 
see [I]). 

Given such an analytic solution, it could then be expanded about various limits. 
An appropriate post-Newtonian expansion would yield the Maclaurin spheroids in 
the zeroth order, then the analytically known first post-Newtonian corrections to the 
next order (see [21 [3]) and so on. Such a post-Newtonian expansion can be expressed 
entirely in terms of elementary functions at least up to the fourth order (4) . Similarly, 
an expansion with respect to the angular velocity would result in the inner and outer 
Schwarzschild solution in the zeroth order. In this paper we shall consider the first 
order contribution to the slow rotation expansion using the formalism introduced by 
Hartle [5]. 

Chandrasekhar and Miller [5] already considered this problem in 1974 by solving 
the equations numerically. Since there now exist computer programs capable of 
solving the complete Einstein equations in the case of stationarity and axial symmetry 
to extremely high accuracy, a further numerical treatment of the slowly rotating 
approximation would only be useful inasmuch as it provides a fairly simple alternative 
for arriving at a good approximation to the full field equations. An analytic treatment, 
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however, can pinpoint where the stumbling blocks preventing further progress lie, and 
can augment the disc solution [T] in suggesting a "lower bound" for the complexity of 
the full-solution. A discussion of some of the analytic issues that arise in this context 
can be found in [7] and references therein. 

We present the basic equations for a slowly rotating homogeneous perfect fluid in 
Sj2] and show that the equation of primary importance is one of Heun's equations. 
In Sj3] we discuss some of the properties of Heun's equations and consider what 
transformations could result in a simpler equation. Various series solutions are 
discussed and derived in SJH and use is made of the numerical solution of the full 
Einstein equations in order to determine an appropriate truncation order. Using these 
truncated series, approximate expressions for all metric functions are discussed in $5] 
and provided in | Appendix A| 



2. Basic Equations for Slow Rotation 



As in L 6] , we devote this section to listing the fundamental equations derived by Hartle 
and then specialize them to the case of constant energy density. 

The metric describing a slowly rotating, stationary and axisymmetric fluid is 
given by 



ds 2 



with 



? 2vs {\ + 2h) dt 2 +e 2Xs 



1 + 



-2m 



dr 2 



+r 2 (l + 2k) dO 2 + sin 2 6 (dtp- u dt) 2 +0{fl 3 ) 



(1) 



h= h (r) + h 2 (r)P 2 {cos9) 
m = m Q (r) + m 2 (r)P 2 (cos9) 
k= k 2 (r)P 2 {cos8) 
uj = uj(r). 

In the above equations, P 2 (cos9) denotes Legendre's polynomial of order 2, O is the 
angular velocity (see [5] for an account of what 'slow rotation' means in terms of f2), 
u> is of order f2 and h, m and k are of order O 2 . In the above expansion in spherical 
harmonics, use was made of a coordinate freedom in order to transform away the term 
Mr). 

The non-rotating metric, obtained when O — > 0, is 

ds 2 = ~e 2us dt 2 + e 2As dr 2 + r 2 (d9 2 + sin 2 6 dtp 2 ) . 

For a non-rotating perfect fluid with pressure ps, constant energy density e, radiusH 
rs and mass 

M S = ^nsrl (2) 
the metric functions ar^l (see e.g. [8]) 
r > rs '■ e' 



r 



r - 2M S 



r < rs 



r - 2M S ' 
1 



B(rY 



-\A-B(r)] 



(3) 
(4) 



2 The coordinate radius of the star rg is not to be confused with Rg in [6], which they use to denote 
the Schwarzschild radius. 

3 We adopt units in which c = G = 1. 
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A := 3W1 



1 



2M s r 2 



(5) 



and the normalized pressure is given by 
Ps = A/3-13 
e B-A ' 

Integrating the equation for hydrostatic equilibrium, it is natural to define 

p* = ln(e+p) - / - d£ 

J e'+p(e') 

= ln(e + p) (e — constant), 
which we expand as with the metric coefficients 

P* = Ps + s Po( r ) + 6 PU r ) ^2(cos6»). 
Hydrostatic equilibrium is then realised when 

I 

3 1 



(0) 



SpQ + ho — —r 2 e 2vs uj 2 = constant = 7 



5p* 2 



1 



(7) 



- 2 -2va~2 n 

-r e b u> = 0, 



where we have introduced^ 

uj := Q — us. 

Letting 

j .— e -(>is+<^s) 
v 2 ~ h 2 + k 2 , 

and restricting ourselves to the case of constant density, we can write the field 
equations as 



d 



du> d 



4<I> d 



dr 2 dr dr ^ ^ r dr 



dmo 
dr 

dh 
dr 

dh 2 
dr 



12 



4 -2 

r J 



dr 



m e 



2A S 



1 



3 dr 



1 



12 



duj\ ' 



— + 8tt Ps - —r 6 f \-r) + 47rr (e + p s ) Sp* 



dr J 



(86) 
(8c) 



+ 



-2e 2As ldv s 
r 2 I dr 

r 2 j 2 ( dw 



V 2 



-2 



du s 1 /dvs ( 1 d(j 2 ) dX s S 



dr 



dr \ 2j 2 dr 



dr 



G 



dr 



dr 



dr dr 



m 2 



re 



,-2As 



1 dv s 
r dr 

-h. + W 

6 



- 4 -2 

-r j 



1 dvs 
dr 

dcu 
dr 



„2A 



s / di>s 
dr 



-r uj — \ — 

3 dr 



duj\' 
~dr~) 



1 3 - 2 l(il 

3 dr 



(8d) 
(8c) 
(8i) 



4 Our notation u> = Q — u> corresponds to lu in [5] and tv in [6]. In this paper, a bar above a symbol 
will denote a dimensionless quantity. 
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In the vaccum region, where e = p = and j = 1, the solutions to the above 
equations are 



r > rg 



2 7 

n - 4 (9a) 



7 2 

m = £M g (96) 

<5M J 2 

^ = 7^2M S + r3(r-2M s ) ( ° C) 

«* = -^-^ , 2MS Qif^-l) (9e) 

m 2 = i r A 2 (-/ l2 + ^) ! (9J) 

where is the associated Legendre function of the second kincjl and J, SM and 
K are constants, J being the angular momentum and 5M the change in mass with 
respect to the non-rotating configuration. 

Equations 1)8 aft are listed in a hierarchical order, each equation being soluble once 
the preceding equations have been solved (note that (|8dj) and (|Sej) form a coupled 
system). Eq. (|8ap plays a pivotal role, since the solution of the remaining equations 
relies on the solution of this equation. Hence we shall devote the next section entirely 
to this equation. Before doing so, it will be convenient to put (|8a[) . which is a Heun 
equation, into the standard form (see e.g. [9]). 

To bring it into this form, let us begin with a transformation of variables in (|8a[) 
from r to the B{r) of ([5|). Making use of the constant A from (J5J), we find 

d 2 uj , . „o . „ n duj 
( _ + (4^-5^ + 1)- 

With the further substitution (see [TU] (2.329)0 
1 - B 



(B - A)(B 2 - 1) — j + (4B 2 - 5AB + 1)— - 4A£ = 0. (10) 



we arrive at 



with 



' : ' 1 (11) 



dz 2 \z 2 — 1 z — a J dz z(z — l)[z — a) 
1 -A 



q = -2A, a = 3, = 0, 7=-, 
5 

<^=2' e = a + /3 + 1 - 7 - 6 e = -1. 

5 The minus sign in front of the second term of ( 19 el ) has to do with the normalization of the Legendre 
function and the choice of the branch cut. 

6 Both this substitution and the resulting (IIH are not unique, since there is freedom as to how the 
original four singularities of the equation, (1, — 1, A, co), are mapped onto (0, 1, a, co). This freedom 
will be considered amongst the transformation of [|3] 
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3. Heun's Equation: Properties and Transformations 

Using the terminology of [9], we consider three categories of solutions to Heun's 
equation: local solutions, Heun functions and Heun polynomials. 

Local solutions to Heun's equation are valid in the neighbourhood of one 
singularity and are associated with one of the two exponents there. Since there are four 
singularities, there are a total of eight such solutions. Heun functions are solutions 
valid in an region containing two adjacent singularities. Heun polynomials are 
solutions to Heun's equation valid at three singularities. Despite the name, not all 
such solutions are polynomials, but are of the simple form 

Hp (z) = (z - If 2 (z ~ a)^p n {z), (12) 

where p n {z) is a polynomial of degree n and 01,2,3 is one of the exponents associated 
with the singularity at z = 0, 1, a respectively. 

In order to determine which type of solution we need to look for, we have to 
determine the domain of relevance for our particular Heun equation given by (jlip . 
The parameter A and the quantities that depend on it are defined on the intervals 

A €(1,3), a € (0,-1), q € (0,-3/2). 

The limit of inifinite central pressure is given by A — > 1 and the Newtonian limit by 
A — ► 3. We exclude these limits in this analysis (as indicated by the open intervals), 
but note that Chandrasekhar and Miller include an interesting discussion of the limit 
of infinite central pressure in [5] and believe that consideration of such limits in the 
context of the confluent Heun equation could lead to interesting results. 

The interior of the star, i.e. the region of validity of (fTTj) . is given by the interval 

a + 1" 



B e [l, A/3] =^ze 



6 



0. 



(13) 



where z = represents the centre of the star and z = (3 — A)/6 its surface. 

Since a is always strictly negative whereas z is always real and positive, we are 
interested in a solution including the singularity at zero and extending along the real 
axis toward the singularity at z = 1 (though not reaching it). The function ui must 
remain finite (and non-zero) at the point z — 0, whence we require at the very least 
the local solution there corresponding to the exponent o~\ = 0. 

The remainder of this section is devoted to a survey of attempts that can be 
made to find a closed form solution to our Heun equation satisfying these constraints. 
Although none of these attempts proved successful, the following subsections are 
intended to stimulate new ideas in this direction, explain where problems lie in what 
may seem promising solution strategies and help the reader avoid choosing fruitless 
paths of thought. 



3.1. Heun Polynomials 

The most tractable solution to Heun's equation is one of the Heun polynomials, so 
that we begin by searching for such a solution even though the validity at the two 
'additional' singularities is not necessary for our purposes. As was mentioned in 
footnote [6j the particular form of the Heun equation given in (|11|) is not unique 
since various transformations map a Heun equation onto a new Heun equation with 
other parameters and exponents. Indeed there are 192 mappings of Heun's equation 
onto itself, corresponding to the 192 local solutions that can be provided for it. 
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Table 1. The values of the exponents cr; of (1 1 2 P and the parameters a and f} 
corresponding to the four classes of Heun polymials that would yield permissible 
solutions to our physical problem. The numbering of the classes follows the 
conventions in [9]. 



Class 




<T2 


&3 


a 




P 




I 











—n 


7 + 5- 


hM 


-n-1 


111 





l-S 





S-n-1 


7" 


he4 


- ra 


V 








1 - e 


e — n — 1 


7" 




- n 


VII 





1-8 


1 - e 


<5 + e - n - 2 


7" 


-n- 


r 1 



These are made up of 24 transformations of the independent variable, which map 
the four singularities onto themselves, and 8 elementary power transformations of the 
dependent variable. A chapter on these transformations is contained in [5] and a 
discussion of their group structure as well as a useful table containing all of them can 
be found in [11] . These transformations will not occupy us further in this paper, but 
the statements made and conclusions drawn are valid (modulo obvious modifications 
resulting from the transformation) for all mappings of onto itself. 

The exponents o~i (see ()12j) ) of the four Heun polynomials consistent with our 
physical problem as well as the parameters a and j3 for the solutions are listed in 
Table [T] A comparison with (fTl) shows us that only a Heun polynomial of class I is 
consistent with the given parameters since n, the degree of the polynomial p n {z) in 
(|12p. can of course only be a non- negative integer. The set of equations that must be 
solved in general in order to determine p n { z ) provides a polynomial equation of degree 
n + 1 for the accessory parameter q. In our case, the equation is simply q = 0, which 
lies just outside the allowed range for q and thus there exists no Heun polynomial, 
which is a solution of the physical problem being considered. 

3.2. Heun Functions 

Since we require a solution valid at z = and extending along the real axis toward 
z = 1 , a Heun Function seems to be the appropriate solution. It is related to the two- 
point connection problem treated in [T2] and discussed further in [T3] . The search for 
such solutions generally implies solving an eigenvalue problem and is valid only for a 
restricted set of values for the accessory parameter q. In our case, q is not independent 
of the singularity a, and no such solutions exist. 

3.3. Rational and other Transformations 

The next step we take in looking for a solution of our Heun equation is to consider 
a broader set of transformations. Since Heun's equations are closely related to 
the hypergeometric equations, one can hope to find a transformation relating the 
two equations. Kuiken [M] addressed the question as to when a hypergeometric 
equation can be transformed into a Heun equation via a rational substitution of 
the independent variable (excluding the trivial case a(3 = q = 0). Her work was 
reexamined and completed by Maier [15]. Using their results, we have found that no 
rational transformation of the independent variable maps our Heun equation (with 
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the specified parameter range) onto a hypergeometric equation^. Even had such a 
transformation been found, it would only have been valid for specific values (or one 
specific value) of the parameter A, because of the nature of the necessary conditions 
for such a transformation to exist. 

Having ruled out rational transformations, we turn our attention briefly to integral 
transformations derived by Carlitz and Valent and described in [16] , which impose no 
restriction on the accessory parameter. Such transformations have proved useful in 
finding previously unknown, closed form solutions to Heun's equation. In our case, 
these methods only allow for marginal progress. Consider first a transformation of the 
independent variable 

p , B-a 
B(l -a)' 

Applying to this the transformation (9) of [16] . and denoting the new quantities by a 
douple prime, one finds that 

7 " = e" = i and a" + 0" = 6" 

and can apply the quadratic transformation discussed on pg. 59 of |17j . Here, however, 
one is not guaranteed that the solution will be regular at the centre of the star. 

More modern ideas, such as relating the Heun equation to the Schrodinger 
equation through the "generalized associated Lame" (GAL) potential seem tailor made 
to the situation being considered here. Indeed, the analogue of (32) in [TS] with b = 3/2 
can be applied to (fTTj) . but the restriction to the accessory parameter is q = as it 
was in the case of a Heun polynomial. 

We know in fact that no meromorphic solution to our Heun equation exists (for 
arbitrary q) since the necessary condition (=l/2-ra,meZ [19] is not met. 



4. Series Solutions 

The power series solution (Frobenius solution) is the most obvious local solution to 
Heun's equation. It could be used to generate a solution regular at the point z = 
and corresponding to the exponent zero, 

oo 

S '.= S c m z , 

m=0 

which converges, in general, within a circle extending out to the next singularity. It 
is also possible to generate a series solution in terms of hypergeometric functions as 
described in [20] 

oo 

H := 2_j c ™^™> witn ( 14 ) 

rn—0 

T(a - 5 + m + 1)T(8 - 6 + m + 1) m , . 

= - p7 , ~n rl l z m F(a + m,0 + m;a + (3-5 + 2m+l;z) : 

I (a + p — o + 2m + 1) 

where T(x) is the gamma function and F(a, b ; c ; z) the hypergeometric function. The 

coefficients c m are given by the three term recurrence relation 

L c + AVi = 

K m c m-i + L m c m + M m c m +i = (m = 1,2,3,.. .) 

7 Many of the computations in this paper made use of the computer algebra program Maple™. 
Maple is a trademark of Waterloo Maple Inc. 
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Table 2. The error of the partial sums Sn and H n is shown at the point on the 
star's surface z = (3 — A)/6 for the value A = 5/2. 



n \S-S n \ \H-H n \ \H- J ^(zH n + H n+1 ) 






2 x icr 1 


2 x icr 1 


2 x 10" 2 


1 


2 x icr 2 


4 x 1CT 3 


2 x icr 4 


2 


2 x icr 3 


i x icr 4 


3 x 1CT 6 


3 


2 x icr 4 


7 x 10~ 6 


6 x 10~ 8 


■1 


2 x icr 5 


5 x 10~ 7 


i x icr 9 


5 


2 x icr 6 


5 x lCi- 8 


5 x 1CT 10 



where 

Krn+l 



(a + m) (/? + m) (e + m) (a + f3 — 5 + m) 



(a + /3-5 + 
= am(7 + m — 1) 



2m)(a + /3-<5 + 2m + l) 
(a + m)(a -<5 + m+l) + (/3 + m)({3 - 



5 + 771 + 1) 



1 



M„ 



a + /3-(5 + 2m-l 
a/3(7 — 2m) — em(S — m — 1) 
a + /3-<5 + 2m + l 
(a — i5 + m)(/5 — <5 + m)m(7 + rn - 



(a + /3 - <5 + 2m - l)(a + - S + 2m - 
— m(a + (3 — 5 + m) — q 



1 



The analysis in [2 j 



2m- l)(a + /3-<5 + 2m)' 
| shows that -ff converges in the domain 



l-(l-«) x /a 



1 - (1 + Z)V2 



< 



l-(l-o) 



1/2 



(16) 



l-(l + a)V2 ■ 

Taking into account the range of z (Eq. [T3|) . in turns out that both of these series will 
converge for ^4 e (3/2, 3), but not in the range A £ (1,3/2]. If one chooses to use this 
approximation as a rough model for neutron stars, one is unlikely to be interested in 
A < 1.5 however 

There are two reasons for opting for the series of hypergeometric functions H 
in favour of the power series S. For one thing, it generally converges more quickly 
than the corresponding power series. In addition, it is an alternating series in which 
successive partial sums straddle the limiting value of the sum, so that they provide 
both an upper and lower bound. A comparison of the two series can be found in 
Table [2] for the value of the function Q at the star's surface z — (3 — A)/6 for A = 5/2. 
The partial sums H n tend to be about one order of magnitude more accurate than 
the corresponding S n . Taking into account that H is an alternating series, a weighted 
average such as the one presented in the last column of the table, ultilizes H n and 
H n+ i to produce a somewhat more accurate estimate for H than would be given by 
H n+ i alone. The weight chosen in the table is based on the fact that each term in the 
series (TlH) is multiplied by z m . 

Whereas the table shows that H converges faster than S for A = 5/2, Fig. Q] 
shows how close and H§ have come to their asymptotic values for the full range 
of convergent A values. We see that both series converge very quickly for A — > 3 and 



8 For a density of e ■ 
2.8M . 



10 18 kg/m 3 , A = 3/2 corresponds in the static model to Mg fa 5.5 X 10 30 kg : 
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10 -io 



10 -ao 

1.5 2.0 2.5 3.0 

A 

Figure 1. \(S-S S )/S\ and \(H-H S )/H\ on the surface of the star z = (3-A)/6 
are depicted over the range A £ (3/2, 3). 

that the hypergeometric series converges faster than the power series for most of the 
convergent interval. 

5. The Approximate Metric 

5.1. Determining the Truncation Order 

Having found a convergent series solution to Heun's equation (for a certain range 
of A), we can now compare the slow rotation approximation with the solution of 
the full Einstein equations for axially symmetric, stationary, homogeneous, uniformly 
rotating fluids to motivate how to choose the truncation order of the partial sum 
H n . The solution of Einstein's equations can be found numerically, and using the 
spectral program described in [21j . we can reach machine accuracy and thus have an 
absolute measure for the accuracy of the Hartle approximation. We choose to compare 
configurations with the same total mas^l M and angular velocity f2, where the bar 
indicates that the quantity is expressed in terms of units of the energy density s. 

We are interested in solutions of Einstein's equations that are connected to the 
non-rotating limit via a continuous parameter transition. Such solutions were studied 
in [22] and termed the 'generalized Schwarzschild' class of solutions. A sequence of 
stars from this class with constant mass M < 4/9-Rs ~ 0.145 can be followed from 
the non-rotating limit to the mass-shedding limit, at which a cusp forms along the 
equatorial rim. Along any such sequence, Q reaches a global maximum somewhere 
between the two limiting points. In the approximation considered here, however, the 
mass is monotonic in A G (3/2,3) for a given f2. Therefore it is not possible to find 
more than one configuration for given M and f2 and we are forced to restrict our 
attention to the portion of the sequence between the static limit and the maximum 
value for fi. 

In Figs[2]and[3]we see how J depends on f2 for one sequence close to the Newtonian 
limit (M = 0.01) and one highly relativistic configuration (M = 0.1). We see verified 

9 A discussion as to how one can integrate 118 al l to determine the mass follows in § 15.21 
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J 





0.0 0.3 0.6 0.9 

n 

Figure 2. The angular momentum J is shown as a function of angular velocity Q 
for a sequence of star's with the constant mass M = 0.01. The solid line depicts 
numerical results, which are accurate to better than 10 digits and thus act as an 
absolute standard of reference. 

0.005 



J 





0.0 0.2 0.4 0.6 0.8 1.0 

Figure 3. As in Fig. [2] but for a sequence with the mass M = 0.1. 

in each plot that the approximation becomes arbitrarily good for O — s- 0. Furthermore, 
we see in agreement with Fig. [T] that H converges quickly for large A, and indeed H§ 
is bairly distiguishable from Hi in this figure (A « 2.8 over the whole range of the 
plot). This does not mean that the slow rotation approximation is good however. The 
discrepancy between the correct and approximated values for J are quite noticable for 
O = 0.6, although its maximal value along this sequence is fi ss 1.23. The inaccuracies 
due to the slow rotation approximation are evident in Fig. [3] as well. Here, however, 
the improvement brought about in going from Hi to H5 is more pronounced (A varies 
from about 2.7 to 2.8 over the range of the plot). 

To have more quantitative information as to when to truncate the series H, 



numerical 
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Table 3. The relative error in the angular momentum as a function of H n is 
shown for various values of Cl for a sequence with M = 0.01 



J-J(Hi) 



J-J(H 2 ) 
J 



J-J(H 3 ) 
J 



0.01 


1.5 x 10" 


2 


2.5 x 10" 


4 


2.6 x 10" 


5 


1.8 x 10 


0.1 


1.7 x 10~ 


2 


1.6 x 10" 


-3 


1.8 x 10- 


3 


1.8 x 10 


0.2 


2.3 x 10" 


2 


7.1 x 10" 


-3 


7.4 x 10" 


3 


7.3 x 10 


0.1 


4.5 x 10- 


2 


2.9 x 10" 


-2 


3.0 x 10" 


2 


3.0 x 10 


0.8 


1.4 x 10" 


1 


1.3 x 10" 


1 


1.3 x 10~ 


1 


1.3 x 10 


1.23 


4.7 x 10" 


1 


4.6 x 10" 


1 


4.6 x 10" 


1 


4.6 x 10 



Table 4. As in Table |3] but for a sequence with M = 0.1. 



J-J(Hi) 




J-J(H 2 ) 




J-.7{H a ) 




J 




J 




J 





J-J(Hoo) 



0.01 

0.1 

0.2 

0.4 

0.8 

1.45 



5.2 x 10~ 2 

5.3 x 10~ 2 
5.8 x 10~ 2 
7.5 x 10~ 2 

1.4 x 10" 1 
4.0 x lO" 1 



1.3 x 10~ 2 
1.1 x 10~ 2 
6.7 x 10" 3 
1.1 x 10~ 2 
7.9 x 10~ 2 
3.6 x 10- 1 



5.8 x 10" 3 
7.3 x 10" 3 
1.2 x 10" 2 

2.9 x 10" 2 
9.5 x 10~ 2 
3.7 x 10" 1 



1.5 x 10" 5 
1.5 x 10~ 3 
6.0 x 10~ 3 
2.4 x 10~ 2 
9.0 x 10" 2 
3.7 x 10- 1 



Tables [3] and 2] consider the same sequences and show the relative error in the angular 
momentum as a function of H n for various values of Q. Each table shows a range of O 
values extending up to its maximum for the respective sequence. The extremely large 
errors emphasize how poor the slow rotation approximation is for homogeneous bodies 
near the mass-shedding limit. One should keep in mind however that these errors are 
not astrophysically relevant since the fastest known pulsars have an angular velocity 
of f2 « 0.4, meaning that the approximation holds to within about 3%. If one chooses 
the order of truncation, by requiring that the relative error in J be in the vicinity 
of a few percent whenever the slow rotation approximation allows for such accuracy, 
then H2 seems an appropriate choice. When providing explicit expressions for metric 
functions in the next section, we shall choose this truncation order. 

5.2. Deriving Explicit Expressions 

If we label the arguments in the hypergeometric functions of H (see (O) by a, 6 
and c, then we find a + b — c = — 3/2. By applying Gauss' relations for contiguous 
hypergeometric functions (see e.g. § 2.8 in [23]), we can convert this into two functions 
with c = —1/2 and c = 1/2. Quadratic transformations are known to exist for such 
functions so that they can be converted into associated Legendre functions of the form 
Pv 71+1 1 2 1 n 6 Z (see e.g. §7.3.1 Eqs (36) and 41 in [24]). Associated Legendre functions 
of this form can, in turn, be represented with a finite number of terms, § 3.6.1 in |23J. 
The explicit expression that results for H2 can be found in the appendix. 

The requirement that the metric functions be continuous at the star's boundary 
leads via (Pa[) to 
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where we have chosen to express the relation in terms of dimensionless quantities. 
Taking into account 

M S = ^Rl <=^Rl = A(i _ A 2 /9), (18) 

(cf. ([2])), we thus know J as a function of Cq and A for a given f2. Since the normal 
derivative of u> must be continous on this surface (see e.g. § 30.5 in [5]), 

d 



dr 



_ 6J 



(19) 



must hold as well. If we were to prescribe A as in [6], then this would suffice to 
determine cq and J. 

If on the other hand, we choose to prescribe ft and M as in the previous section, 
then we have to solve (|86f before being able to determine both A and c$. 

Eq. (l86l) can then be integrated from the centre out to the surface of the star to 
give an expression for SM = m | r= ^ s + J 2 /Rg from J9&J. The constant of integration 
must be chosen such that mo(0) = so that g rr remains finite at the star's centre. 
The total mass of the system is then M = M s + 5M = 4/3?ri?| + SM. Taking into 
account our choice of normalization (cf. @) 

M S = ^nRl ^R 2 s = J^(l - A 2 /9), 

we have arrived at expressions containing the unknown variables c and A, which can 
then be determined by prescribing M and Cl. 

Because the integrand in the integral to determine mo contains quadratic terms 
in lu and its derivative, it does not seem possible to obtain an analytic expression for 
mo based on the hypergeometric representation discussed above. In the appendix, 
expressions are thus provided for the metric functions in terms of a series in z. 

Eq. (|Fc[) for determining ho is simply a first order, linear differential equation, 
where ([7]) is used to replace Sp^ by an expression in hy. The constant of intergation 
can be chosen arbitraril\F°l and is chosen in the appendix such that /io(0) = 7. This 
choice amounts to identifying a rotating body with a non-rotating one of the same 
central pressure since the choice implies through (|7|) that 6po(0) — 0. The nature of 
such choices of identification and some of the effects different choices can have was 
discussed in the context of the post-Newtonian approximation in [J. The interior 
solution for ho must join continuously to its exterior solution, thus fixing the constant 
7 of 0. 

The system of equations (|8rf[) and (|8ej) for obtaining hi and V2 can be solved as 
described in [6] . Both functions must vanish at the origin so that h and k have unique 
values there. The remaining constant as well as the constant K in (|9<ij) and (|9ej) can 
then be determined by requiring that these two functions be continuous at the star's 
boundary. 
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through the SFB/TR7 "Gravitationswellenastronomie" and by a "Landesgraduierten- 
stipendium" from Thuringia. 



Appendix A. Explicit Representation of the Truncated Metric 



The formulas for the metric functions are collected in this appendix. They are valid 
over the range A £ (3/2, 3) and make use of the variable 



k = 



The symbol F(a, b; c; z) refers to the hypergeometric function. The metric depends 
on three paramters (one scaling parameter and two "physical" parameters) that were 
here chosen to be A — 3yT — 2Ms/rs (see ©), rs and Co- Ms and rs are the mass 
and radius (in Schwarzschild coordinates) of a spherical star and Co is an integration 
constant (see (|14|) and lfTB"|) ) giving the value of oj at the star's centre w(0) = 4y/ir/3co. 
The expressions below do not diverge as one approaches the Newtonian limit A — > 3 
since cq — > sufficiently quickly in that limit. 



4^7? c 



1 



8A; 



b(A-l 

which can be evaluated to yield 



8{Az) 



10 = 



((A-i)zy(z-i) 



13 



-1/24- 



23 



35(^1-1) 2 
z - 63) A 2 



M 2,5; y;z 



Zs/z 



2/3 



V~zA 



arcsin 



(Vi) 



21 



z A z -\ z A 



45 45 ' 12 
- 2/9 z (4 z 2 - 10 z + 3) A + 4/3 z 2 (z - 1) 
or expanding as a power series and dropping higher terms 



CO 



4^/Trco 



1 



8A 



&A(9A - 



5(A-1) 35(A-1) 2 



We now leave off the ellipsis at the end of each expression. 



m 



5l2Trc 2 r s 3 Az 5 / 2 



840 (A - l) 2 + 12 (A - 1) (A - 125) z 



SM = 



525 (9 - A 2 f /2 {A- I) 

+ (2975 + 602 A - 441 A 2 ) z 2 

16Wr s 3 (3-A) 5/2 
127575 (9-A 2 ) 3/2 [A- if 

V6A (84015 - 109128 A + 34922 A 2 + 3176 A 3 - 441 A 4 



+-| (9 - A 2 ) (A - 1) (A + 3) 5/2 (45 - 56 A + 9 A 2 ) 2 ) 

Do / 



h = 7 - 



1792Avr c 2 r s 2 
15(yl-l) 3 (A 2 - 9) 



(1.1a) 



(1.16) 



(1.1c) 



(1.2) 



(1.3) 
(1.4) 
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with 

-9SM 647rc 2 r s 2 (A-3) 



7 



A 2 r s ' 893025A 2 (A + 3) 01 - l) 4 

x (81 A 8 - 522 A 7 - 2102 A 6 + 14262 A b + 65961 A 4 
-165591 A 3 - 47466 A 2 + 298890 A - 164025) 



, f (3A-25)/ 512Wrg^4 , 

ft2= /OZ+ 3T-J TT o Z 

I 7(A - 1) 3 (A 2 -9) (A- if 1 



2 



' 2(5A 2 -17A+ 104)/ 5127rc 2 rs 2 ^ (129 A- 575) \ 
21 (A - l) 2 225 (A 2 - 9) (A — l) 4 i 



z 3 (1.5) 



-2/ 5127rc 2 r s 2 A ^ ^ 2 

4(A- 13)/ 10247TCQ 2 r s 2 A (7/4 + 25)" 



y 2 = -; a z 

V^-l 3 (A 2 - 9) (A - 1) 



7(A-1) 2 75 (/4 2 - 9) (A - l) 4 

-2(5 A 2 -26 A + 221) / 



z 3 



\ 21 (A -if 

5127TC0 2 r s 2 A (51 A 2 - 2262 A + 9275) 
1575 (A 2 - 9) [A - if 



z 4 (1.6) 



where the continuity of /i 2 and V2 at the star's boundary then give for the constant /o 
and K of E3 and (El 



/o = irc 2 Rs s x 



-512 (/4 - 1) (A 2 - 9) (81 A 8 - 522 A 7 - 3560 A 6 + 23658 A 5 



I 2 (H. 



f- 44370 A 4 - 120771 A 3 - 2046924 A 2 + 3774195 A - 492075) Q2 (^ ~ 1 

- 768 A (l62 A 9 - 1206 A 8 - 3160 y4 7 + 32728 A 6 + 14001 A 5 - 328536 A 4 

f- 373590 A 3 - 2968596 /4 2 + 4333095 A + 32805o)q 

28350 (A - 3) (A + 3) 2 (A - if 

x (1989 - 2514 A + 962 A 2 - 74 A 3 + 5 A 4 ) (j^- - 1 j 

- 14175 (A - 3) (A + 3) (A - l) 2 A 

(5 A 4 - 92 A 3 + 1790 A 2 - 5052 A + 4149) Q\ f -0- - 1 



(1.7) 
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K 



7rco 2 rs 3 A 



64 (2835 A 13 - 65331 A 12 



893264 A 1 



8037225r s (A — 1) (A + 3) 
4933304 ,4 10 - 8773582 ,4 9 + 172596150 .4 s - 361180413 A 7 - 958944555 A 6 
1549090818 A 5 + 11997038790 A 4 - 41313070791 A 3 + 54020240631 A 2 

33095488275 A + 8082331875 



2 (A — 1)(A+ 3) (5 A 4 - 74 A 3 + 962 A 2 - 2514 A + 1 



^ (5 A 4 - 92 ,4 3 + 1790 A 2 - 5052 A + 4149) Q\ [ - 1 



M s 



Ms 



m,2 can then be calculated algebraically using |8j 
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